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Sondaae a ponderation spatiale d'une voie de transmission 
La presente invention concerne un procede de sondage 
d'une voie de transmission formee de plusieurs canaux. II 
s'agit done d'une technique connue sous le nom de diversite 
5 de-reception , technique selen laquelle un r6cepteur comporte 
plusieurs antennes chacune associee h un canal de 
transmission distinct. En d'autres termes, 1* invention 
propose une methode d f estimation de la reponse 
impulsionnelle de ces canaux de transmission. 

10 En effet, dans un systeme de transmission, notamment 

par ondes radio, le r6cepteur regoit pour chaque canal de 
transmission un signal emis par un emetteur. Le signal emis 
subit des fluctuations d' amplitude et de phase dans le canal 
de transmission, si bien que le signal regu par le r6cepteur 

15 ne lui est pas identique. Les fluctuations du signal sont 
essentiellement dues & ce que I'homme de metier appelle 
1 'interference intersymbole. Cette interference peut 
provenir de la loi de modulation employee pour la 
transmission et elle est egalement due k la propagation 

20 multi-trajets dans le canal. 

II se trouve que le signal regu est g6n6ralement issu 
d'un grand nombre de reflexions dans le canal, les 
diff6rents trajets empruntes par le signal emis conduisant 
ainsi k des retards varies au niveau du recepteur* La 

25 reponse impulsionnelle du canal repr6sente 1' ensemble de ces 
fluctuations, auxquelles est sounds le signal emis. II 
s'agit done Ik de la caracteristique fondamentale 
representant les transmissions entre 1' emetteur et le 
r6cepteur . 

30 La reponse impulsionnelle du canal est utilis<§e 

notamment par un egaliseur qui a pr6cisement pour fonction 
de corriger 1 1 interference intersymbole dans le recepteuf. 
Une methode classique pour realiser une estimation de cette 
reponse impulsionnelle consiste & disposer dans le signal 

35 emis une sequence d • apprentissage formee de symbol es connus. 
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Cette sequence est choisie en fonction cie la loi de 
modulation et de la dispersion du canal, dispersion devant 
s* entendre ici comme le retard d'un symbole emis erapruntant 
le trajet le plus long du canal par rapport a ce meme 
5 symbole empruntant le trajet le plus court. La dispersion 
est couramment exprim6e comme un multiple de la duree qui 
separe deux symboles 6mis successifs, soit un nombre de 
"duree symbole 11 . 

A titre d'exemple, on citera deux techniques connues 
10 d' estimation de la r6ponse impulsionnelle d'un canal de 
transmission . 

La premifere technique fait appel a des sequences 
d 1 apprentissage particuliferes dites sequences CAZAC, pour 
l 1 expression anglo-saxonne "Constant Amplitude Zero 

15 Autocorrelation". De telles sequences sont d6crites dans 
1* article de A. MILEWSKI : "Periodic sequences with optimal 
properties for channel estimation and fast start-up 
equalization"/ IBM Journal of Research and Development, 
vol.27, N # 5, Sept. 83, pages 426-431. 

20 Le systfeme de radiocommunication cellulaire num6rique 

GSM fait appel & des sequences d 'apprentissage TS form6es de 
26 symboles not6s a c h &25 prenant la valeur +1 ou -1. Ces 
sequences poss&dent les proprietes suivantes : 

20 

25 £ a ± 2 = 16 

i = 5 

20 

2 a i a i + k = 0 si 0 < |k| £ 5 

i = 5 

En no tan t d la dispersion du canal qui vaut 4 dans le 
cas du GSM, 1 1 estimation de la reponse impulsionnelle prerid 
30 la forme d'un vecteur X & 5 composants notes Xq & x 4 . 

La sequence de symboles reqrus S correspondant a la 
sequence d 1 apprentissage TS est form§e elle aussi de 26 
symboles not6s s Q a S25. On suppose naturellement ici que 
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1 1 emetteur et le recepteur sont parf aitement synchronises et 

dans ce cas l 1 estimation de la reponse impulsionnelle X est 

donnee par 1' expression suivante : 

20 

xfc = — X) ai s l + k pour 0 £ k s 4 

16 i = 5 

5 La technique CAZAC pr§sente l'avantage d'une grande 

simplicity de raise en oeuvre. Cependant, on remarque que 
chaque composante de la r6ponse impulsionnelle est 6tablie a 
partir de seulement 16 symboles re^us. Etant donn§ que la 
sequence d • apprentissage comprend 26 symboles et que la 
10 dispersion du canal vaut 4, il y a de 1 ■ information dans le 
signal regu qui n'est pas prise en compte et cela conduit a 
line reduction des performances par rapport a l'id€al 
th6orique . 

La deuxi&ne technique connue fait appel au crit&re des 
15 moindres carr€s. Elle est d6crite notamroent dans les 
demandes de brevet FR 2 696 604 et EP 0 564 849. En matifere 
de rappel, cette technique fait appel & une matrice de 
mesure A construite & partir de la sequence d 1 apprentissage 
TS de longueur n. Cette matrice comprend (n-d) lignes et 
20 (d+i) colonnes, d representant toujours la dispersion du 
canal. L* 6 lament figurant & la ifeme ligne et a la jieme 
colonne est le (d+i-j)i£me symbole de la sequence 
d 9 apprentissage : 



25 



f &4 a3 a2 ai ao 
as a4 a3 a2 &i 
a6 as a4 a3 a2 
A = 37 

*25 a 21 
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La sequence d'apprentissage est choisie telle que la 
matrice A fc A soit inversible oil l'operateur .t represente la 
transposition. C'est naturellement le cas pour les sequences 
CAZAC mais c'est egalement le cas pour d'autres sequences, 
5 Dans la sequence de symboles rejus, on ne prend pas en 

compte les quatre premiers sq & S3 car ceux-ci dependent 
egalement de symbole inconnus §mis avant la sequence 
d'apprentissage , 6tant donn§ que la dispersion du canal vaut 
4. Par un abus de langage on d§finira done dorSnavant le 

10 signal regu comme un vecteur S ayant pour composantes les 
symbol es regus, S4, S5, sq, ... , S25. 

Des lors, l 1 estimation de la r6ponse impulsionnelle 
prend la forme suivante : 

X = (A fc A)" 1 At . S 

15 Cette technique des moindres carr£s est un peu plus 

complexe que la pricSdente mais il faut noter que la matfice 
(At A)" 1 At est calcul§e une seule fois. On remarque ici que 
chacune des composantes de 1' estimation de la r§ponse 
impulsionnelle X est etablie a partir de 22 symboles re£us 

20 et non pas de 16 comme dans le cas de la technique CAZAC On 
doit done s'attendre & une amelioration des performances. 

Cependant, quelle que soit la technique utilis6e, les 
reponses impulsionnelles de chacun des canaux de la voie de 
transmission sont consider6es comme independantes . 

25 La pr£sente invention a ainsi pour premier objet un 

proc6d6 de sondage d f une voie de transmission qui prend en 
compte le fait que les diff6rentes antennes sont 
spatialement li6es. 

Selon l 1 invention, le proc€d§ de sondage d'une voie de 

30 transmission s 1 applique lorsque cette voie est form§e de 
plusieurs canaux et il n§cessite une estimation de la 
reponse impulsionnelle Ci, C2, . C n de ces canaux. Le 
proc6d<§ comprend les stapes suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
35 transmission, 
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- etablissement d'une reponse impulsionnelle corrigee au 
moins par ponderation de ces estimations de reponses 
impulsionnelles au moyen de cette statistique spatiale et 
d'une estimation du bruit additif de ces canaux. 
5 Avantageuseraent, cette statistique spatiale correspond 

a une estimation de la correlation des canaux de 
transmission pris deux & deux. 

Suivant un mode de realisation pr6f erentiel , cette 
estimation de la correlation des canaux de transmission 

10 prend la forme d'une matrice de correlation spatiale ainsi 
definie que l 1 element figurant & la ifeme ligne et a la jifcme 
colonne soit obtenu par lissage du produit de la transposee 
hermitienne de la reponse impulsionnelle estim6e du ieme 
canal et de la reponse impulsionnelle estim6e du jieme 

15 canal . 

Selon une caracteristique additionnelle, un signal 
regu S par un premier canal correspondant a une sequence 
d'apprentissage emise, en notant A la matrice de mesure 
associee a cette sequence d'apprentissage, 1' estimation du 
20 bruit additif Nqi de ce canal est obtenue par normalisation 
de 1 1 energie du vecteur ( S - ACj ) . 

On peut pr6voir que cette normalisation soit suivie 
d'une etape de moyennage. 

De plus, une matrice de bruit N etant form6e & partir 
25 de 1' estimation du bruit additif Nqi/ Nq2# •••/ N 0n des 
canaux, une matrice de ponderation spatiale 6 9 etant ainsi 
definie a partir de la matrice de correlation spatiale G et 
de la matrice de bruit : 

G 1 = G (G + N)" 1 , 
30 les reponses impulsionnelles corrigees C'i, C f 2/ ••• / C' n 
sont obtenues par 1' expression suivante : 
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Le proc6d6 de sondage d'une voie de transmission est 
done mis en oeuvre & partir de 1 1 estimation de la reponse 
impulsionnelle des diff6rents canaux considers comme autant 
5 de canaux independants . Or les erreurs d' estimation sont 
inevitables. De manifere g6nerale, la determination de la 
r6ponse impulsionnelle d'un canal de transmission unique est 
un problfeme qui ne peut §tre r6solu de fagon exact e en 
presence de bruit additif. De plus les techniques 

10 ant6rieures font l'hypothfese implicite que cette reponse 
Impulsionnelle peut prendre une forme quelconque. 

La pr6sente invention a ainsi pour deuxifeme objet une 
methode d' estimation de la reponse impulsionnelle d'un canal 
de transmission qui presente une meilleure resistance au 

15 bruit additif ou, autrement dit, qui mfene a une erreur 
r edui te compar 6e & 1 1 erreur d 1 estimation des techniques 
connues. Cette methode sera avantageusement appliqu^e pour 
l»un au moins des canaux formant la voie de transmission 
avant 1'etablissement des reponses impulsionnelles corrig6es 

20 par ponderation au moyen de la statistique spatiale et d'une 
estimation du bruit additif de ces canaux. 

Selon l 1 invention, cette methode necessite un signal 
regu par un canal, ce signal regu correspondant a une 
sequence d'apprentissage 6mise. La methode comprend les 

25 etapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission, 

- etablissement de 1' estimation de la r6ponse impulsionnelle 
de ce canal , cette estimation 6tant ponderee par la 

30 statistique temporelle du canal au moyen du signal regu. 
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La statistique temporelle du canal repr£sente une 
valeur de la reponse impulsionnelle ant6rieure & 
1 'acquisition du signal re^u. La susdite pond6ration 
introduit le fait que la reponse impulsionnelle aff^rente au 
5 signal legu a une valeur probabiement pltrs piroche de cette 
valeur ant6rieure qu'une valeur qui en serait trfes 61oign6e. 
Ainsi, statistiquement , l'erreur d' estimation est diminu6e. 

Avantageusement , cette statistique correspond a une 
estimation de la covariance de cette r6ponse impulsionnelle. 
10 Selon une premiere variante du proc6d6, celui-ci 

comprend les stapes suivantes : 

- lissage de la reponse impulsionnelle et orthonormalisation 
au moyen d'une matrice de transformation W pour obtenir 
l f estimation de la covariance qui prend alors la forme d , une 

15 matrice L f , 

- recherche des vecteurs propres v^ 1 et valeurs propres Xj/ 
associ6es de cette matrice L 1 , 

- - estimation de la reponse impulsi onnel 1 ft instantanee du 
canal a partir du signal regu et application de cette 

20 matrice de transformation W pour former un vecteur X», 

1 1 6tablissement de I 1 estimation pond6ree Xp etant ainsi 
realisee : ^ 



= Z 



r ii_^£( v ; h .x.)l 



1 h 
W V ± n 



oil Nq est un nombre r6el strictement positif repr§sentant le 
25 bruit additif du canal. 

On peut ici pr^voir que le bruit additif soit choisi 
6gal & la plus petite des valeurs propres Xj/. 

Lorsqu'un sous-ensemble de ces valeurs propres Xj/ 
prSsente une contribution inf6rieure & un seuil 
30 pr6d6termin6 , on peut Sgaleraent pr6voir que. chacune, de ces 
valeurs propres soit forc6e a la valeur du bruit additif. 

La complexity s'en trouve r6duite d*autant. 

Selon une seconde variante du proc6d§, 1 9 estimation de 
la covariance prenant la forme d f une matrice R, en notant A 
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la matrice de mesure associee a la sequence d'apprentissage, 
1 'etabl is semen t de 1 1 estimation ponder^e est ainsi realist : 

Xp = (A*- A + NqR" 1 )" 1 At.S 
oCi No est un reel strictement positif representant le bruit 
5 additif du canal. - . 

Par ailleurs, il est possible d'effectuer une etape 
d'orthonormalisation de la matrice R au rooyen d'une matrice 
de transformation W pour obtenir une nouvelle matrice R 1 , 
1* estimation pond#r£e prenant alors la nouvelle forme 
10 suivante : 

Xp » W fc (I + NqR 1 " 1 )" 1 wt A't.S 
oh la matrice A 1 est egale au produit de la matrice de 
transformation W et de la matrice de mesure A. 

Avantageusement, l 1 expression (I + NqR 1 " 1 )" 1 est 
15 calcul€e au moyen du lemme d 1 inversion matricielle. 

La pr6sente invention apparaitra maintenant de manifere 
plus d6taillee dans le cadre de la description qui suit oil 
sont proposes des exemples de mise en oeuvre & titre 
illustratif, ceci en reference aux figures annex6es qui 
20 repr6sentent : 

la figure 1, un diagramme identifiant les 
principales etapes d'un mode de mise en oeuvre du procede 
de sondage d'une voie de transmission selon I 1 invention , 

- la figure 2, une premiere variante d 1 estimation de 
25 la r6ponse impulsionnelle d'un canal selon l 1 invention, 

- la figure 3, une deuxieme variante d 1 estimation de 
la r6ponse impulsionnelle d'un canal selon l 1 invention, 

L* invention sera pr6sent6e dans son application au GSM 
car ce systfeme a le mSrite d*§tre bien connu de I'homme du 
30 metier. II s'agit done Ih d'une presentation adoptee dans un 
souci de clart6 mais il ne faut y voir en aucun cas une 
limitation de l 1 invention & ce seul syst&ne. 

En r6f6rence h la figure 1, le procide de sondage 
d'une voie de transmission s 1 applique lorsque cette voie 
35 comporte au moins deux canaux de transmission, n canaux dans 
le cas g6n6ral. Les canaux correspondent chacun a une 
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antenne distincte. On considere done que chaque canal a fait 
1'objet d'une estimation de sa reponse impulsxonnelle 

C n selon l'une quelconque des 



10 



15 



• • / 



respective C lf C 2 , 
techniques disponibles . 

En-premier lieu le procede prevoit 1 • acquisition d'une 
statistique spatiale de cette voie de transmission. Par 
statistique spatiale on entend un ensemble de donnees 
refletant le comportement de cette voie sur une perxode 
predetermine que l'on conviendra d'appeler periode de 
correlation. En effet, les differentes antennes etant fxxes, 
les signaux recus sur ces antennes presentent un certaxn 
degre de correlation. L' invention entend precisement tirer 
parti de ce fait pour ameliorer la qualite d' estimation de 
la reponse impulsionnelle de l'un des canaux au moins. A 
titre d'exemple, une telle statistique peut etre obtenue au 
moyen d'une matrice de correlation spatiale G : 



rv ( c h Cl ) V(C h C 2 ) V( C hc n ) 1 



G = 



VCC^Cj) v(c h c 2 ) 



v( c h c n ) 



[v{c* Cl ) v( C hc 2 ) 



v( c{;c n ) 



20 



25 



ou l'operateur . h repr6sente la transposition hermitienne. 

Cette matrice G carree de dimension (n,n) peut ainsi 
etre representee de maniere generique par 1* element gi j 
figurant a la ieme ligne et a la jieme colonne : 

«ij = V(C i C D ) h 
cet element gij est obtenu par lissage du produxt C i C j 

au moyen" des reponses impulsionnelles estimees C if Cj des 
ieme et jieme canaux obtenues durant la periode de 
correlation. Ce lissage est une estimation de la correlation 
de ces deux canaux. 
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On entend ici lissage dans un sens tres general/ 
c'est-a-dire toute operation permettant de lisser ou de 
moyenner le produit C^C^ sur la periode de correlation. 

Un premier exemple de lissage consiste a faire la 
5 moyenne de ce produit sur cette p£riode suppos^e comprendre 
e estimations successives de chacune des reponses 
irapulsionnelles Cj, C2, ♦ ../ C n : 




1 



10 Un deuxifeme exemple de lissage consiste a actualiser, 

a la pifeme estimation regue pour chacun des ifeme et jieme 
canaux, l 1 expression de lissage V p _ ^C^Q. ) obtenue a la 

(p-l)ifeme estimation au moyen d f un coefficient 
multiplicatif a, ce coefficient compris entire 0 et 1 6tant 
15 souvent connu sous le nom de facteur d'oubli de lissage : 

v p< c i c j > - ac i c j + d-a)V p _ j^CcJcj ) 

L 1 initialisation peut se faire par tous moyens, 
notamraent au moyen de la premifere estimation obtenue ou bien 
par une moyenne obtenue comrae dans le premier exemple pour 

20 les premieres estimations regues. 

Le proc6d§ de sondage propose ensuite de rSaliser une 
estimation du bruit additif N 0 i, N 02/ . .., N 0n present dans 
chacun des canaux au moyen des estimations des r6ponses 
impulsionnelles respectives C lf C 2 f c n de ces canaux. 

25 II sera propose diff6rentes solutions pour mener 6 

bien cette estimation du bruit et on traitera le cas d»un 
seul canal, le premier par exemple, 6tant donn6 que les 
mfemes solutions s'appliquent a chacun des canaux. 

Une solution simple consiste & affecter N 0 i d'une 

30 valeur pr6d6termin6e qui reflate un seuil en dessous duquel 
il est peu probable que le bruit additif puisse descendre. 
Cette valeur pourrait §tre d6terminee par une raesure de 
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rapport signal S bruit, ou par les performances du 
recepteur, ceci a titre d'exemple. 

Par ailleurs, 1' estimation du bruit N 01 peut etre 
obtenue au moyen de l 1 estimation de la reponse 
5 impulsionneUe Cx du premier canal et du signal regu 
correspondant S. Cette estimation du bruit a done pu §tre 
r6alisee avant la mise en oeuvre du procede de sondage h 
ponderation spatiale selon quelque technique que ce soit. 
Toutefois, si ce n*6tait pas le cas, on propose ici une 
10 m6thode appropri6e lorsque 1' estimation de la reponse 
impulsionnelle est acquise selon la technique des 

moindres carr6s. 

II s'agit done d'estimer directement le bruit additif 
a partir du signal regu S et de la matrice de mesure A. En 
15 effet, en notant Ni le vecteur bruit affectant le signal 
regu, il vient que : 

S = AC X + N x 

Compte tenu du fait que les vecteur s S et Nj ont 22 
composantes, l 1 estimation du bruit Nqi prend la forme 
20 suivante : 

Nqi = (^) (S - AC!)h (S - ACi) 

Naturellement cette estimation du bruit additif Nqi 
peut §tre moyenn6e ou liss§e. 

A partir des estimations du bruit additif des 
25 diff6rents canaux N 0 i/ Nq2/ N 0n* on construit une 

matrice de bruit N dont l 1 Element figurant a la ieme ligne 
et a la jifeme colonne est not6 n^ j . 

Lorsque i est 6gal & j, cet 616ment nij prend la 
valeur du bruit estim6 Noi* 
30 Lorsque i et j sont diff6rents, l'§16ment nij est nul. 

Par- - ailleurs^- tou jours dan& un— souci de 
simplification, il est possible de calculer la valeur 
moyenne N 0 des bruits estimSs N 0 i/ Nq2/ — * N 0n et de 
forcer chacun des 616ments nij de la diagonale de cette 
35 matrice N a cette valeur moyenne Nq. En notant I la matrice 
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identite, il s'ensuit que la matrice de bruit N prend la 
forme suivante : 

N = N 0 I 

La matrice de correlation G et la matrice de bruit N 
5 permettent de d6finir une nouvelle matrice, la matrice de 
ponderation spatiale G 1 : 

G« = G(G + N)" 1 
Les reponses impulsionnelles corrig§es C'i, C*2/ •••/ 
C* n au moyen de cette matrice de ponderation spatiale G 1 
10 sont ainsi d6finies : 







f Clt ) 
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oil l f op6rateur repr§sente tou jours la transposition. 

L'une au moins des reponses impulsionnelles corrig€es 

15 C'i est utilisee h. la place de l f estimation de la rgponse 
impulsionnelle dans le r6cepteur. 

A titre optionnel, 1' invention permet 6galement 
d'amSliorer 1' estimation de la r6ponse impulsionnelle C^, 
C2# C n de chacun des canaux, ceci avant la mise en 

20 oeuvre du proc6d6 de sondage a ponderation spatiale d'une 
voie de transmission, c'est-&-dire avant 1 f 6tablissement de 
la ou des r6ponses impulsionnelles corrigSes C 1 1 , C ' 2 / •••/ 

On prSvoit ainsi 1 'acquisition d'une statistique 
25 temporelle pour l'un au moins de ces canaux, le premier par 
exemple. Par statistique temporelle, on entend un ensemble 
de donnSes reflStant le comport ement de ce canal consid§r6 
independamment des autres canaux sur une p6riode 
pr§d6termin6e que l'on conviendra d'appeler p^riode 
30 d 1 analyse ♦ II s'agit done d'une representation du 
comportement moyen du canal pendant la p€riode d 1 analyse. 
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Cette statistique peut etre etablie par quelque moyen que ce 
soit et en quelque lieu que ce soit. En effet, 
1 •etablissement de la statistique peut prendre place dans un 
equipement quelconque du reseau de radioconununication . Ce 
5 -qui iinporte c'est que le r£cepteur puisse acqu^rir cette 
statistique . 

A titre d'exemple, une telle statistique temporelle 
peut §tre obtenue de la manifere suivante. 

Au moyen du signal regu S pendant la p6riode 
10 d f analyse, on calcule une estimation X de la r6ponse 
impulsionnelle selon une m6thode connue. 

Si I'on retient la technique des moindres carr6s cette 
estimation X vaut : 

X = (A* A)-l At.S (1) 
15 On rappelle ici que l*6metteur et le r6cepteur sont 

supposes synchronises h mieux que un demi-symbole pres, 
auquel cas le signal de r6ception est form6 par le vecteur S 
dont les composantes sont les symbol es regus S4 h S25 
synchrones des symbol es a4 a a 2 s de la sequence 
20 d'apprentissage TS. si une telle synchronisation n'§tait pas 
acquise, plusieurs solutions sont disponibles pour 
l'acquSrir et on en citera deux a titre d'exemple. — 

La premiere solution consiste & dScaler le signal regu 
en avance ou en retard de j p6riodes symboles de sorte que 
25 Sj fc = (s 4 _j, s 5 _j, s 6 -j, ... s 2 5-j). 

On calcule alors 1 1 estimation Xj pour chacun des 
vecteurs Sj et on retient la valeur jjj pour laquelle Xj^*Xj 
est maximal. Cette valeur donne la synchronisation 

escompt€e et il suffit de remplacer le vecteur S dans 
30 I 1 equation (1) par le vecteur Sjjj. 

La deuxi&me solution consiste & augmenter 
artificiellement la dispersion d du canal d'une quantity 2q 
pr6determin6e . On peut alors d€finir une matrice de mesure 
modifi#e A m comprenant (n-d-2q) lignes et (d+2q+l) colonnes. 
35 En reprenant les valeur s 26 et 4 respectivement pour n 
et d : 
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Am = 



' a 4 + 2q 
a 5 + 2q 
a 6 + 2q 
a 7 + 2q 



a 4 a 3 a 2 a l a 0 
a 5 a 4 a 3 a 2 a l 
a 6 a 5 a 4 a 3 a 2 



a 25 



a 21 - 2q 



II faut 

alors r6duire le nombre de composantes du signal regn S de 
cette raeme quantity 2q et on peut convenir de retenir le 
5 vecteur S m modif i6 : 

Sin* = ( s '4+q/ s's+q,... s' 2 5-q) 
On obtient ainsi vine estimation modifiee X m : 

*m = (Am* A m ) Ajnt-SjQ 
Cette estimation modif i6e X m comporte d+2q+l 
1 0 composantes : 

Xj^^ = (x_q/ Xq / / •••/ X4 f • •• X4+q) 

On recherche alors la valeur jjj de j comprise entre -q 
et +q qui maximise 1' expression suivante : 



k = 0 

15 oii l'op€rateur .* repr6sente la conjugaison complexe. 

Cette valeur jjj determine 1 • estimation X de la r6ponse 
impulsionnelle pour une dispersion d=4 : 

xt = (Xj M , Xj M+1/ ... Xj M+4 ) 

La synchronisation s'en d6duit imm6diatement en 
20 appliquant le dScalage j M au signal regu S. 

On peut dfes lors appliquer & nouveau !• Equation (1). 
On construit maintenant une matrice de lissage L par 
lissage des differentes estimations X obtenues pendant la 
pSriode d 1 analyse pour obtenir une estimation de la 
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covariance associee & cette reponse impulsionnelle. On 
entend ici lissage dans un sens tres general, c'est-&-dire 
toute operation permettant de lisser ou de moyenner la 
r§ponse impulsionnelle sur la periode d' analyse. On obtient 
5 ainsi , una representation statistique du comportement du 
canal de transmission, Ce lissage peut etre obtenu notamment 
par l'une des deux m6thodes proposees ci-dessus, 
1' expression la plus simple de cette matrice de lissage L 
etant la suivante : 
10 L = ^ ^ XX h 

ou m correspond au noxnbre de sequences d'apprentissage sur 
lequel est calculi le lissage. 

On admet ici que cette matrice de lissage peut £tre 
approch6e par liquation suivante : 

15 L - (AtA)" 1 N 0 + R (2) 

oil Nq represente tou jours le bruit present dans le 
canal de transmission ou bruit additif et oil R est une 
matrice que I f on a coutume d'appeler statistique a priori du 
canal car elle represente le comportement du canal 

20 abstraction faite du bruit. 

On admet egalement que la matrice de mesure A est bien 
conditionnee , c'est-^-dire que les valeurs propres cle la 
matrice A fc A sont tr£s proches les unes des autres . Dans ce 
cas, il est int6ressant de proc6der & 1 1 orthonormalisation 

25 des vecteurs constitu6s par les lignes de la matrice de 
mesure A, mais il ne faut pas voir 14 une limitation de 
1* invention. 

Pour ce faire, on emploie une matrice de 
transformation W telle que : 
30 A 1 = AW et A^A' = I 

oil I represente la matrice identity. 

En notant L' la matrice ainsi definie : 
L = WL'Ytft, 

il vient que liquation (2) peut maintenant s'6crire : 
35 L 1 ~ N 0 I + R 1 (3) 
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Selon une premifere variante representee a la figure 2, 
on remarque que les vecteurs propres Vi 1 de L 1 et v^ de R' 



^ de R' sont decal6es de Nq. Soit en prenant tou jours la 
5 meme valeur de 4 pour la dispersion du canal, pour tout i 
compris entre 0 et 4 : 



10 propres et valeurs propres de R' et celle de L 1 sont 
identiques sous reserve que Nq soit connu. 

L 1 etape d 1 estimation du bruit sera dScrite plus loin 
pour rendre 1* expose plus clair, bien que cette etape 
precede celle qui va maintenant §tre explicitee. 

15 il convient done maintenant de rechercher les couples 

valeur-propre/vecteur-propre pour l'une ou 1* autre des 
matrices L 1 ou R 1 . Cette etape ne sera pas plus detailiee 
car bien connue de l'horame du metier. Par ailleurs, il va 
sans dire que l'on peut annuler les valeurs propres dont la 

20 contribution est jug6e non significative • Par exemple, si 
ces valeurs propres sont class6es par ordre decroissant, on 
supprime les dernieres qui sont t el les que leur somme soit 
inf6rieure & un seuil predetermine. 



25 impulsionnelle instantanee X etablie selon l'une quelconque 
des techniques connues & partir du signal regu correspondant 
a la dernifere sequence d'apprentissage regue. En notant 
X = WX 1 , cette dernifere estimation est pond6r6e par la 
methode suivante pour obtenir une pond6ration temporelle Xp 

30 de la reponse impulsionnelle instantanee : 



sont identiques tandis que les valeurs propres X^ 1 de L 1 et 



Xi* = X ± + N 0 

II apparalt ainsi que la determination des vecteurs 



L 1 etape suivante consiste a estimer la reponse 
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Pour obtenir la pond6ration Xp il faut done estimer le 
bruit additif Nq. 

Cette estimation du bruit peut §tre effectuee 
notamment selon I'une des m§thodes exposees plus haut. 
5 Une autre possibility consiste & consid6rer que la 

dernifere valeur propre, (la plus faible) de la matrice de 
lissage L est 6gale & Nq : 

A, 4 » = N 0 ou X4 = 0. 
Quelle que soit la m6thode retenue, la pond£ration 
10 teraporelle Xp de 1 1 estimation de la rSponse impulsionnelle 
instantan6e peut alors etre r#alis6e comme mentionn6 ci- 
dessus • 

Selon une seconde variante representee h la figure 3, 

on etablit !• estimation pond6r§e Xp directement comme suit : 

15 Xp = (A fc A + NqR" 1 )"^^ 

ou bien en reprenant la matrice de transformation W 

d6firiie ci-dessus : 

^ Xp = W (I + NqR'" 1 )" 1 W fc A't.S (4) 

Conforro6ment a 1" Equation (3) : 

20 R 1 = L 1 - N 0 I 

II convient Ih aussi d 1 estimer le bruit additif N 0 . 

Une solution avantageuse pour obtenir la pond6ration 

temporelle Xp consiste & adopter la m6thode suivante. 

On divise la matrice R' par Nq : 
R • 

25 B = 

N 0 

II s'ensuit que : 

I + NqR 1 " 1 = I + B" 1 
On utilise le lemme d 1 inversion matricielle pour 
calculer la matrice de pond6ration P = (I + B" 1 )" 1 . 
30 Ainsi, en notant e^ les vecteurs canoniques, on 

proc&de h 3J iteration-suivante % 

- initialisation : 

P = B 

- pour i variant de 0 h d (4 dans le cas present) : 
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P 6tant connu, il reste a etablir la pondSration Xp 
selon 1' equation (4). 

On remarquera que la matrice de pond6ration P n'est 
5 pas necessairement calcul6e a 1' occasion de chaque nouvelle 
sequence d'apprentissage 6mise. Elle peut Stre calculee a un 
rythme plus lent car elle varie sensiblement au meme rythme 
que R 1 et done plus lenteraent que le signal re$u S. 

On remarquera 6galement que 1 1 estimation ponder§e est 
10 r6alis6e sans avoir recours a la r6ponse impulsionnelle 
instantanee. Elle est produite directement k partir du 
signal regu S. 

Quelle que soit la variante retenue, e'est done 
1» estimation pondir6e Xp qui sera avantageusement utilis6e 
15 coinme estimation de la reponse impulsionnelle Cj pour la 
mise en oeuvre du proc6d6 de sondage a pondSration spatiale 
d'une voie de transmission , e'est-a-dire pour produire une 
ou plusieurs r§ponses impulsionnelles corrig6es. 
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REVEND I CAT I ONS 

1) Procede de sondage d'une voie de transmission 
formee de plusieurs canaux, ce procede necessitant une 
estimation de la reponse impulsionnelle (C^ C2, #C n ) de 

5 ces canaux, caracterise en ce qu'il coinprend les Stapes 
suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
transmission , 

- etablissement d'une r6ponse impulsionnelle corrig6e (C^, 
10 C r 2/ • •* * c9 n) au rooins par pondSration desdites estimations 

de r6ponses impulsionnelles (Cj, C2/ ... /C n ) au moyen de 
ladite statistique spatiale et d'une estimation du bruit 
additif (N 01 , N 0 2/ N 0n ) desdits canaux, 

2) Proc6d6 selon la revendication 1, caracterise en ce 
15 que ladite statistique spatiale correspond & une estimation 

de la correlation desdits canaux de transmission pris deux & 
deux. 

3) Proc6de selon la revendication 2, caracterise en 
ce que cette estimation de la correlation des canaux de 

20 transmission prend la forme d'une matrice de correlation 
spatiale (G) ainsi definie que l 1 element (gi j ) figurant a la 
ieme ligne et & la ji&me colonne soit obtenu par lissage du 
produit (C^Cj) de la transposee hermitienne de la r6ponse 

impulsionnelle estimee (Ci) du i£me canal et de la reponse 
25 impulsionnelle estimee (Cj) du ji&me canal. 

4) Procede selon la revendication 3 caracterise en 
ce que, un signal regu S par un canal correspondant & une 
sequence d 1 apprentissage 6mise, en notant A la matrice de 
mesure associee & ladite sequence d 1 apprentissage , 

30 l 1 estimation du bruit additif (Nqi) de ce canal est obtenue 
par normalisation de l*6nergie du vecteur t$ - AC^) . 

5) Procede selon la revendication 4, caracterise en ce 
que ladite normalisation est suivie d'une etape de 
moyennage . 

35 6) Procede selon l'une quelconque des 

revendications 3 a 5 caracterise en ce que, une matrice de 
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15 



20 



25 



30 



bruit (N) 6tant form6e & partir de 1 1 estimation du bruit 
additif (N 01 , Nq2> •••/ ^On) des canaux, une matrice de 
ponderation spatiale (G 1 ) etant ainsi d§finie h partir de 
ladite matrice de correlation spatiale (G) et de ladite 
matrice de bruit : 

G' = G (G + N)" 1 , 
lesdites r£ponses impulsionnelles corrig§es (C 1 !, C*2, , 
C' n ) sont obtenues par 1' expression suivante : 
/c, l t ) (c£ ^ 

= G 



v c 'n J l, c n' J 

7) Proc6de selon l'une guelcongue des 
revendications 1 a 6 caract6ris6 en ce que, avant 
1 1 6tablissement de ladite rSponse impulsionnelle corrigee 
(C*!) d»un canal , le signal regu (S) par ce canal 
correspondant & une sequence d'apprentissage 6mise, il 
comprend les Stapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission , 

- 6tablissement de 1* estimation de la reponse impulsionnelle 
(Cx) dudit canal, cette estimation (Xp) 6tant pond^ree par 
la dite statistique temporelle du canal au moyen dudit 
signal regu (S). 

8) Proc6d6 selon la revendication 7, caract6ris6 en 
ce que ladite statistique temporelle correspond a une 
estimation de la covariance de ladite rSponse 
Impulsionnelle . 

9) Proc6d6 selon la revendication 8 caract6rise en 
ce qu'il comprend les 6tapes suivantes : 

lissage de ladite r€ponse impulsionnelle et 
orthonormalisation au moyen d'une matrice de transformation 
W pour obtenir ladite estimation de la covariance qui prend 
alors la forme d'une matrice L 1 , 
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- recherche des vecteurs propres (v^ 1 ) et valeurs propres 
(A-i 1 ) associees de cette matrice L f , 

- estimation de la rfiponse impulsionnelle instantanee du 
canal a partir dudit signal regu (S) et application de cette 
matrice de transformation W pour former un vecteur X* , 
l'^tablissement de ladite estimation pond§r6e (Xp) 6tant 
ainsi r6alis6e : 

f a • _ » ^ 

h 



= £^LpK h -*.))wv:> 



oil Nq est un nombre reel strictement positif repr6sentant le 
10 bruit additif dudit canal. 

10) Proc6d6 selon la revendication 9, caract6ris6 en 
ce que ledit bruit additif (No) est choisi 6gal & la plus 
petite desdites valeurs propres (X^ 1 ). 

11) Proc6d§ selon I'une quelconque des revendications 9 
15 ou 10/ caract6ris6 en ce qu'un sous-ensemble desdites 

valeurs propres (Xj/) pr6sentant une contribution inffirieure 
a un seuil pr§d6termin§ r chacune de ces valeurs propres est 
forc6e a la valeur dudit bruit additif (No). 

12) Proc§d€ selon la revendication 8 caract6ris6 en 
20 ce que ladite estimation de la covariance prenant la forme 

d'une matrice R, en notant A la matrice de mesure associee & 
ladite sequence d 1 apprentissage , 1 • 6tablissement de ladite 
estimation pond£r€e (Xp) est ainsi r€alis€ : 
Xp = (A fc A + NqR" 1 )" 1 A fc .S 
25 od Nq est un r6el strictement positif reprfesentant le bruit 
additif dudit canal. 

13) Proc6d§ selon la revendication 12, caract6ris6 
en ce qu'il comprend une 6tape d'orthonormalisation de 
ladite matrice R au moyen d'une matrice de transformation W 

30 pour obtenir une nouvelle matrice R* £ cette estimation 
pond6r£e prenant alors la nouvelle forme suivante : 

Xp = Wt(i + NqR 1 " 1 )* 1 W fc A ft .S 
oil la matrice A 1 est 6gale au produit de cette matrice de 
transformation W et de ladite matrice de mesure A. 
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14) Proc6de selon la revendication 13, caracteris6 
en ce que l 1 expression (I+NqR 1 " 1 )" 1 est calcul6e au moyen du 
lemme d 1 inversion matricielle. 
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c 2 



G = 



VCQ^) V(Q h C 2 ) VCCfo) 

V(C 2 h Q) V(C 2 h C 2 ) V(C 2 h Cn) 



^(0*0) 



No 



N (estimation) 



G' = G (G+N)" 1 



CV 



IC'„V 



= G' 



FIGURE 1 
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L'(XX h ) 



L , (XX n ) = N 0 l + R , 



Nq (estimation) 



(v'j.A'j)? i = 0,...,d 
Lv'j = X\ v'j 



X = (A t A)~ 1 A t .S 



WX'= X 



Xp=L^ 2 (v , | 1 X , )WVi 



i=0 /t| 



FIGURE 2 
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L'(XX h ) 



Nq (estimation) 



P = (l + NqR'" 1 )~ 1 



A'= AW 



X n = WPA ,t S 



FIGURE 3 
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